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ABSTRACT 


We  present  here  an  extension  of  Hagstrum's  aodel  for  the  interpretation 
of  ion  neutralization  spectra.  We  take  into  account  the  fact  that  the 
neutralization  process  depends  on  both  the  density  of  states  and  the  orbital 
size  and  takes  place  continuously  at  all  surface-ion  separation  distances  at 
which  the  ion  and  metal  orbitals  overlap.  The  model  is  applied  to  the  calcula 
tion  of  the  density  of  states  of  the  Nl(lll)  surface. 
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I.  INTRODUCTION 

It  has  been  demonstrated  by  Hagstrum  that  ion  neutralization  can  be 

1-4 

used  to  probe  the  outermost  layer  of  a  solid  surface.  This  is  experimentally 
achieved  by  exposing  the  surface  to  a  low  kinetic  energy  (10-100  eV)  ion  beam. 
The  ions  are  neutralized  by  a  two  electron  Auger  process,  which  causes 
electron  ejection  from  the  surface.  The  kinetic  energy  of  the  ejected 
electrons  contains  information  about  the  electronic  states  of  the  surface. 

The  development  of  metastable  quenching  spectroscopy^  ^  has  increased 
the  interest  in  ion  neutralization  spectroscopy  by  extending  it  to  very  low 
(i.e.  ~0.05  eV)  incident  kinetic  energies.  This  technique  exposes  the  surface 
to  a  thermal  beam  containing  ground  state  and  metastable  noble  gas  atoms. 

Under  certain  conditions  (i.e.  when  the  work  function  of  the  surface  is  larger 

or  only  slightly  lower  than  the  ionization  potential  of  the  metastable)  the 
impinging  metascable  is  ionized  by  resonant  charge  transfer  to  the  surface. 

The  resulting  noble  gas  ion  collides  with  the  surface  and  generates  an  ion 
neutralization  spectrum. 

The  Interpretation  of  IN  spectra  is  based  on  a  simple  model  proposed 
2  4 

by  Hagstrum,  *  which  expresses  the  spectrum  in  terms  of  a  convolution  of  the 
electronic  density  of  states  with  itself.  In  certain  cases  it  is  necessary 
to  extend  this  model  to  include  the  matrix  elements  of  the  Auger  neutralization 
process  and  the  motion  of  the  incident  ion.  In  this  paper  we  provide  such  an 
extension,  which  is  simple  enough  to  be  a  practical  tool  for  the  interpretation 
of  the  IN  spectra.  The  model  could  also  be  useful  in  the  analysis  of  other 
surface  science  techniques  in  which  ion  neutralization  plays  a  part :  electron 

g 

stimulated  desorption  (ESD) ,  electron  stimulated  desorption  ion  angular 
distribution  (ESDIAD)»^  secondary  ion  mass  spectroscopy  (SIMS)^*^  and  others. 
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II.  THE  MODEL  DESCRIBING  THE  AUGER  NEUTRALIZATION  PROCESS. 

II.  1.  The  Curve  Crossing  Description. 

The  Auger  process  is  simply  described  by  the  curve  crossing  model  explained 

below.  The  orbitals  involved  in  the  Auger  transition  are  shown  in  Figure  1: 

ii.  is  the  singly  occupied  orbital  of  the  impinging  ion,  tU  and  V,  are  the 

orbitals  containing  the  two  metal  electrons  involved  in  the  transition,  and 

1/^  is  the  metal  orbital  (located  above  the  vacuum)  co  which  the  ej.ected 

electron  is  promoted.  In  the  initial  state  (Figure  la)  there  is  a  hole  in 

'll.  and  iii,  is  emoty.  In  the  final  state  v  and  'it  have  a  hole  and  and  'j, 
i  k  •  U  V  i  k 

are  doubly  and  singly  populated,  respectively 

To  describe  the  energetics  of  this  system  we  use  the  energy  of  the  ground 
state  metal  plus  that  of  the  neutral  atom  at  infinite  distance  from  the 


surface,  as  a  zero  energy  reference.  Within  the  one  electron  picture 
discussed  above  the  initial  state  energy  is  £j^(z) ,  and  the  final  state  energy 
is  -y  V(z).  Here  V(z)  is  the  energy  of  interaction  between  the 

neutral  and  the  ionized  surface  when  they  are  separated  by  a  distance  z, 


and  is  the  energy  of  the  orbital  . 

The  energy  conservation  condition  for  an  Auger  process  taking  place  when  the 
ion  is  at  a  distance  z  from  the  surface  is 


£,  +  e  +  e  +  V(z)  =  £.(z) 
k  u  V  i'  ' 


We  note  that  it  is  not  necessary  to  use  the  one-electron  energies  to  describe 
this  system.  The  quantity  £^(z)  -  V(z)  =  I(z)  can  be  replaced  by  the 
ionization  energy  of  the  atom  located  at  the  distance  z  from  the  surface; 


_ £,  +  £^ + £^  can  be  replaced  by  the  total  energy  of  the  excited,  positively 

+  * 

charged  metal  surface  (denoted  (Me  ).  ,) created  by  the  .Auger  neutralization. 

yvK 

The  dependence  of  the  initial  and  final  state  energies  on  the  particle 


surface  distance  is  schematically  plotted  in  Figure  2.  The  upper  cuf/e , 


called  in  what  follows  the  ionic  cuf/e,  represents  the  interaction  energy 
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of  the  metal-ion  system  (Me+A  ) ,  which  is  dominated  by  the  image  force  at 

longer  distances  and  is  repulsive  at  shorter  ones.  The  lower  curve,  called 

in  what  follows  the  neutral  curve,  is  the  energy  of  the  neutral  atom  inter- 
+  * 

acting  with  (Me  The  neutral  curve  has  a  much  weaker  distance  dependence 

than  the  ionic  one. 

The  separation  between  the  curves,  called  here  the  energy  mismatch,  is 

-  I(2)-£^-a  .  (2) 

We  introduce  the  notation  ct  =  since  the  energy  mismatch  depends  on  a 

only,  not  on  and  separately. 

+  * 

If  we  assume  chat  the  interaction  energy  between  the  atom  and  (Me 
is  independent  of  snd  (at  fixed  £^) ,  then  the  neutral  curve  is  highly 
degenerate:  there  are, 

states  having  different  £,^  and  £^,  hut  Che  same  a.  One  of  the  models  proposed 
4 

by  Hagscrum  assumes  that  Che  Auger  process  cakes  place  at  fix  a’^ao  and  the 
number  of  Auger  electrons  having  energy  £,  is  proportional  to  N  .  . 
Deconvolution  of  the  Auger  spectrum  can  therefore  provide  the  density  of 
states  0(c). 

If  we  consider  two  values  of  a,  namely  a'  and  a"  (and  fix  £^) ,  we  generate 
two  neutral  curves  as  shown  in  Figure  3.  These  curves  have  different 
asymptotic  energy  mi.smatches  A,  .  i  (®)  and  A,  .,(”)  equal  to  I(“)— and 
I(«o)-a''-£^,  respectively.  They  cross  the  ionic  curve  at  two  different  points 
Z(a')  and  Z(a'')  ,  which  can  be  obtained  by  solving  the  equation 

A,^  (Z(ci))  -  I(Z(a))-£, -a=0  (3) 

for  and  i»a". 


Since  Che  Auger  neutralization  process , generating  an  Auger  electron  with 


kinetic  energy  £  can  Cake  place  for  all  the  neiicral  curVes  chat  cross  the 
k  ’ 

ionic  one  (i.e.  for  which  the  equation  ^(z)=0  has  a  solution)  we  must, 
in  principle,  consider  the  whole  family  of  neutral  curves  labelled  by 
Thus  the  ionic  curve  interacts  with  a  continuous  band  of  neutral  curves, 
which  are  schematically  drawn  in  Figure  4. 

II.  2.  The  ionization  race. 


The  simplest  model  for  Che  computation  of  the  race  of  production  of  an 
electron  with  energy  uses  the  Golden  Rule  formula 


W^(a)  =  (2Tr/-ft)  /  d£^  /  d£^,  LH(£^^,£^,£^;Z(a)):^  o(£ Jc(£^;o(£^) 


6(1  -£  -£  -£,  )  6(a-J  -£  ) 
a  y  \j  k  U  v 


(4) 


We  consider  here  only  one  neutral  curve  ot,  which  crosses  the  ionic  cur'/e 

at  z=Z(a),  where  the  ionization  potential  is  1^=1 (Z(a)),  One  of  the 

6-functions  in  Equation  (4)  keeps  only  chose  values  of  £^  and  £^  chat  add 
up  to  a;  Che  ocher  enforces  Che  energy  conservation  specified  by  Equation  (1). 

The  matrix  element  H(£^,£i^,£^^;Z(a) )  ,  for  which  we  use  the  notation  H  when 
an  explicit  designation  of  the  variables  is  not  necessary,  will  be  discussed 
in  Section  II. 3. 

As  the  ion  approaches  the  surface  on  the  crajector/  z=z(c),  it  passes  through 
various  crossing  points  Z(a) ,  corresponding  to  various  values  of  a.  3v  using 

z(c)=Z(a),  where  Z(a)  is  the  crossing  point  for  cur-ze  a,  we  can  find  the  function 

a(c)  giving  the  curve  being  crossed  at  time  c.  The  probability  ?.  (c)  that  an 
electron  of  kinetic  energy  £,^  is  produced  at  time  C  (i.e.  at  Che  position 
z(c)=‘Z(a)  where  the  cur'ze  a(c)  is  crossed)  is  given  by 


?,^Ct)  =  p(z(t))(P^^  Wj^(ct(c))  ; 


(5) 


where  p(r(c))  is  the  probability  that  the  ion  reached  the  point  z(c)  without 
being  neutralized , is  the  probability  chat  the  Auger  electron  escapes  into 


the  vacuum,  and  W^(a(.t))  ia  the  rate  of  electron  production  given  by  Equation  (4). 
The  survival  probability  p(z(t)  is  given  by  the  rate  equation 

dpz(t') 

- jp-  *  -W(a(t))  pz(t))  (6) 

where 

W(a)  »  /  d£^W^(a)  .  ‘  (7) 

Assuming  that  we  know  the  trajectot7  of  the  ion  we  can  write  the  finite 
difference  version  of  Equation  (6)  as 

P(Zi+l)  -  P(z^)  =  -W(a^)p(z^)v‘^  ^^i+l“^i^ 

where  z^  =  z(t.),  a  =  ci(t^) ,  v^  =  dz/dt]^^^  ,  and  the  points  t  form  a  fine 

i  1 

grid  on  the  time  axis. 

The  above  equations  provide  a  computational  scheme  for  the  probability 
of  generating  an  electron  with  kinetic  energy  £^,  after  the  incident  particle 
has  completed  its  trajectory.  To  start,  we  must  provide  an  equation  of  motion 

md^z/dt^  ”  F(z)  (9) 

to  propagate  the  ion  and  generate  its  trajectory.  The  force  F(z)  is  discussed 
in  Section  II. 4, 

To  compute  W(a^),  which  is  needed  in  Equation  (8),  we  must  make  a  model 
for  H(£,  ,£  ,£  ;z  ),  use  it  in  Equation  (4)  to  compute  W  (a.),  which  is  then 

tC  y  V  1  iC  1 

used  in  Equation  (7)  to  compute  W(a^).  We  can  now  solve  Equation  (8)  to  obtain 
the  survival  probability  p(t^),  then  use  an  empirical  formula  for  the  escape 
probability  (P,  and  calculate,  from  Equation  (5),  the  probability  ?.  (t.)  that 
an  electron  of  energy  £^  is  generated  at  t^.  Adding  up  over  all  the  times, 

N.  =  :  p(t  )  ,  (10) 

K  ^  1 

we  obtain  the  probability  that  an  ion  generates  an  electron  of  kinetic  energy 
.  This  is  the  computed  IN  spectrum,  which  can  be  compared  to  the 


experimental  results  (after  removing  from  the  measured  spectrum  the  contribution 
form  the  secondary  electron  emission)  to  determine  either  the  density  of 


states  that  best  fits  the  data  (if  one  has  a  reasonable  model  for  the  matrix 

element  H)  or  to  get  the  matrix  element  (if  the  density  of  states  is  known) . 

II.  3.  The  relationship  to  previous  work. 

Since  the  present  work  is  an  extension  of  that  of  Hagstrum  it  is  useful 

to  show  explicitly  the  connection  to  his  model.  To  do  this  we  must  make 

several  approximations.  First  we  assume  that  the  rate  of  electron  production 

is  dominated  by  one  of  the  curves  of  the  family  labelled  by  namely  chat 

corresponding  to  This  reduces  the  problem  to  a  single  curve  crossing 

(Figure  2) .  This  approximation  is  reasonable  only  if  is  very  small  for 

all  ct<a^  and  has  a  large  increase  at  ,  such  that  the  probability  of  having 

been  neutralized  for  z>Z(ct^) ,  and  the  probability  having  survived  as  an  ion 

for  z<Z(a  ) ,  are  both  minimal, 
o 

Under  these  conditions  the  rate  of  production  of  an  electron  of  energy 
is  proportional  cc 

/  de^  |h!^  0(1^  ““ 

o 

where  Ct  is  Che  only  value  of  a  which  we  need  Co  consider, 
o 

If  we  assume  chac  p(£,  )  is  independenc  of  c,  and  chac  H(€,  ,£  ,£  ;2''a  )) 

iC  X  X  U  V  o 


is  indeoendent  of  £  ,£  and  C,  we  can  write  that  is  proportional  to 
P  V  k  k 

-Ji 


o:  /  dS  o(Ej^+3)c(El.^-3) 


(12) 


with 


£.  =(I(Z(a  ))  -£,  )/2 
ok 


(13) 


Thus  Che  IN  spectrum  generated  by  this  model  is  the  convolution  of  the  density 


of  states  with  itself;  the  use  of  a  deconvolution  program  can  supply  t(C), 


.  -  .3 

LT 


H '  “  and  d(£.  )  deoend  stronzlv  on  c  ,£  and  we  can  dei 


is  measured.  If 
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1/2 


0  (£)  5  lH|  oCe)  pCe) 
k 


(14) 


and  use 


«  /  p  CE^+6)  p(E^-6)d6 


(15) 


Clearly  Equation  (15)  is  less  useful  chan  Equation  (12)  since  it  provides  only 
Che  product  of  various  quantities  which  are  of  interest  individually. 

II.  4.  The  model  for  Che  matrix  elements. 


The  assumption  made  in  Hagscrum  model  chat  the  matrix  element  H  is 


independent  on  and  is  likely  to  he  inadequate,  the  nature  of  the  surface 


orbitals  changes  with  Che  orbital  energy.  Consider,  for  example,  the  case 
of  a  metal  that  has  two  bands  A  and  B:  A  is  narrow  and  its  energy  is  larger 
than  chat  of  broad  band  B;  furthermore,  the  orbitals  of  A  are  localized 
tightly  around  Che  metal  atoms,  while  chose  of  B  reach  further  in  space.  (This 
is  consistent  with  the  fact  that  the  band  A  is  narrower  than  B) .  Such  a 
situation  appears  in  the  case  of  Ni,  for  example,  where  the  3d  band  corresponds 
to  Band  A  and  Che  s-p  band  Co  B,  and  the  3d  orbital  are  much  tighter  chan  the 
s  ones. 

If  we  were  to  apply  Hagscrum  model  to  this  system  we  need  to  consider  two 
distances,  and  Z^, corresponding  to  the  distances  at  which  the  band  A  or  3, 

respectively,  neutralizes  most  efficiently  the  hole  in  the  incident  ion.  Since 
the  orbitals  of  band  B  reach  further  in  vacuum  than  chose  of  A,  the  distance 


Z  is  larger  chan  Z .As  the  ion  approaches  the  surface  it  will  first  reach 
D  A 


Z  and  Che  Auger  neutralization  can  occur  bv  electron  transfer  from  B  to  the  ion 

O 


and  electron  ejection  from  either  3  or  A.  We  denote  these  processes  as  BB  and 
3A  ,  respectively.  The  surviving  ions  can  then  reach  and  be  Auger  neutralized 
by  electron  transfer  from  .A  to  the  ion  and  electron  ejection  from  either  A  or  3 


1 


(processess  .AA  and  .AB) . 


-a- 


The  shape  of  Che  IN  speccrun  is  influenced  by  the  difference  in  Che  spatial 
properties  of  the  orbitals.  At  the  distance  Z  Che  BB  process  generates 

D 

electrons  in  Che  kinetic  energy  range  (See  Figure  5) . 

-  2=2 

and  .AB  in  Che  range 


KZg)  -  £1-2  1  =5,  1  I(Zg)-h-^F  • 

-At  the  distance  the  kinetic  energy  of  the  ejected  electrons  is  in  Che  range 

I(Za)  -  2-2  <  <  I(Z^)  -  2-^ 


for  AA,  and 


"V-=2- 


U2^)- 


c  — c 

"F  "2 


for  BA.  These  ranges  are  schematically  shown  in  Figure  5. 

If,  for  example,  the  density  of  states  of  the  B  band  is  very  high  (and 
much  higher  than  that  of  A)  we  expect  EB  to  be  very  efficient  and  .A3  moderately 
efficient.  Furthermore,  most  of  Che  ions  are  neutralized  at  Z  and  very  few 

D 

reach  Z^.  For  chose  that  do,  neutralization  by  the  BA  process  is  more  efficient 
Chan  -AA.  The  resulting  IN  spectrum  having  a  very  intense  BB  emission,  a  weak 
AB  one,  and  an  even  weaker  BA  one,  is  schematically  displayed  in  Figure  7  as 
a  full  line.  If,  for  example,  the  band  A  reaches  far  in  space  and  has  a  high 
density  of  states,  Che  spectrum  is  dominated  by  the  .AA  process,  has  a  weak 
BA  component  and  an  even  weaker  .A3.  The  resulting  spectrum  is  shown  schemacical 
by  Che  dotted  lines  in  Figure  7.  Clearly  the  shape  of  the  spectrum  is  affected 
by  both  the  density  of  states  and  the  spatial  extent  of  the  orbitals  forming 
the  bands.  The  latter  element  is  missing  in  the  Hagstrum  model. 

■A  quantitative  evaluation  of  the  .Auger  matrix  elements  is  very  difficult 
since  the  relevant  electronic  wave  functions  are  those  of  the  surface  in  the 
oresence  of  the  neutral  or  the  ion.  For  this  reason  we  intend  to  pursue  the 
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Integration  of  Eq.  (A17)  gives  the  following  formula  for  the  IKS 
matrix  element  for  an  incident  ion  at  a  separation  distance  z 
from  the  metallic  surface  plane 
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Now  the  focus  of  attention  shifts  to  the  inte^’ratior  over 
the  variable  p.  We  define  the  following  variable  t  r  an  s  f  o  r -na  t  i  o  n 


.  -2  .  , ,  2,1/2 

t=p  andg=(t*u) 


A 1  4  a  ) 


and  thus 


dt  =  2pdp  and  dg  = 


2  1/2 

dt/[2(t'u‘')^^‘^] 


A 1 4b ; 


which  gives 

dt  =  2gdg 


A14c  ) 


Thus  Che  integral  over  the  variable  p  becomes 


ca  '•22.1/2  0® 
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Xdppe^'^  =Xdgge^” 


(  a:3  ) 


Integration  of  Eq.  (A15)  by  "parts"  gives 
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Substituting  the  result  of  Eq.  (A16)  into  Eq.  (A13),  we  obtai.o 


the  following  integral  e.xpression  for  the  INS  matrix  element 


We  begin  to  evaluate  Eq  .  (A9)  by  concentrating  on  the 

integral  over  the  variable  p.  We  first  define  the  following 
transformation  of  variables 
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and  therefore 
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dq  =  p ( p  *5 )  dp  . 
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Through  the  use  of  Eq.  (AlO),  the  integral  over  the  variable  p 
becomes 


J*  dq  e 


-Xq  _  e 


-X( z-z ' -u ) 


for  X  >  0 


We  proceed  by  evaluating  the  integral  over  the  variable  z 
incorporating  the  result  given  in  Eq .  (All) 
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Substituting  Eqs.  (AlO)  through  (A12)  into  Eq.  (A9),  we  obtain 
the  following  e. repression  for  the  INS  matri.x  element 


metal  surface  completes  the  integration  over  the  variable  r  ^  and 

gives  the  screened  Coulomb  interaction  between  the  metallic 

conduction  band  wave  function  il)  (r.)  and  the  ground  state  ion 

pi'  1 

wave  function  t  ^  )  at  the  point  Q  with  the  metallic  conduction 
band  wave  function  over  the  entire  vacuum  region 
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We  now  proceed  with  the  integration  over  the  variable  r^. 
We  must  contend  with  both  the  spherical  geometry  of  the  ion 
ground  state  wave  function  and  the  "planar"  geometry  of  the 
metallic  conduction  band  wave  functions.  Once  again  the 
cylindrical  coordinate  system  is  the  system  of  choice.  Thus  the 
form  of  the  integration  over  the  variable  r^  is  identical  to  the 
previously  employed  in  the  integration  over  the  variable  r^-  Th 
integral  expression  for  the  I. VS  matrix  element  thus  becomes 
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2,1/2 


where  y  is  introduced  to  truncate  the  integral  as  the  separation 


distance  between  the  ion  core  and  the  solid  surface  approaches 


1  3 


width  dp  located  on  this  z  =  z'  surface.  The  surface  element 
contained  within  this  annulus  is  {2np)dp.  All  points  residing 
within  this  annular  surface  element  are  equi-distant  from  the 
axial  point  Q.  Thus  the  screened  Coulomb  Interaction  between  the 
component  of  the  metallic  electronic  wave  function 

contained  within  the  annulus  and  the  elements  of  the  ground  state 
ion  wave  function  ^  f  )  and  the  metallic  electronic  wave 
function  at  the  point  Q  is  given  by 
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Wliere  C  is  defined  C.Cj^C^C^.  Since  we  approx  -nate  the  wave 
function  of  the  excited  metallic  electron  to  be  constant  over  all 
space,  its  presence  in  the  INS  matrix  element  formula  is  denoted 
simply  by  Cj^.  The  screened  Coulomb  interaction  between  the  wav'e 
functions  and  point  Q  and  the  metallic  wave 

function  contained  in  the  entire  plane  z  =  z'  is  obtained 

by  integrating  Eq.  (AS)  over  all  the  annuli  of  which  the  surface 
ix.y.z')  is  composed 
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Fig.  11.  The  wave  function  for  the  incident  helium  ion  is  an 
exponentially  decaying  function  maximized  at  the  position  of  the 
ionic  nucleus 

'Ij,  =  exp(-3r)  .  (.A3) 

The  electronic  wave  functions  of  the  metallic  conduction  band 
have  been  modeled  to  decay  exponentially  in  the  vacuum  region 
according  to  the  following  form  =  Cj  expl-cr^z)  for  z  >  0. 

Thus  points  of  equivalent  metallic  wave  function  amplitude  form 
planes  parallel  to  the  surface  of  the  sera  i  -  i  n f  i  n  i  t e  metal  slab. 

We  focus  on  the  plane  located  at  z  =  z' .  The  amplitude  of  a 
metallic  wave  function  at  any  point  in  the  (x.y.z  =  z')  plane  is 
denoted  iJPj  =  Cj  exp(-ajZ').  The  geometry  associated  with  the 
metallic  wave  functions  suggests  the  use  of  a  cylindrical 
coordinate  system.  In  this  cylindrical  coordinate  system,  the 
wave  function  for  the  Incident  helium  ion  becomes 

il),  =  C.  exp[-p(p^  -  (A4) 

while  form  of  the  metallic  wave  functions  remain  unchanged.  Thus 
the  screened  Coulomb  potential  acquires  the  form 

-*-»  ?  ?1'2  7  712 

Wlr^-r^)  “  exp{ -X[p"-(  z-z  ' -(J)  "]  /  [  p  "  -  (  z  -  z  ' -u  )  "  ]  ' ' 

(  A5  ) 

We  consider  an  infinitesimally  thin  annulus  of  radius  o  and 
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APPENDIX  A: 


An  Outline  of  the  Derivation  of  the-  INS  Matrix  Eleaent 


In  order  to  evaluate  the  integrals  appearing  in  Eq  .  (131 

(and  thus  the  matrix  element  associated  with  Auger 
neutralization),  we  consider  a  physical  system  composed  of  a 
semi-infinite  metallic  slab  occupying  the  half-space  (x.y,z<0) 
and  an  incident  noble  gas  ion  at  the  position  (0.0, z). 

The  matrix  element  eljaracterizing  the  Auger  neutralization 
process  is 

=  /  d"?^  /  dr^  'll  ^  i  r’j  ) (  r  ,  )  W  { "r  ,  -  ir  ^  )  >1'^  (  > 'I’y  '  ' 

(All 

where  integration  is  performed  over  the  soacial  coordinates  of 
both  the  neutralizing  electron  and  the  ejected  electron 

("r^l-  interaction  coupling  the  two  participating  electrons 

is  modeled  by  a  screened  Coulomb  potential  of  the  Yukawa  form 

Wlr"^  "^2*  '  exp(-X  "r.-Tj  ) -' (  ^i"^2  '  ' ' 

where  the  inverse  screening  length  associated  with  the  metallic 
electron  density  in  the  vacuum  region  is  denoted  by  the  symbol  X 
A  schematic  representation  of  the  physical  system  is  displayed  i 
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1 

quenching  spectroscopy  can  generate  with  ease  IN  spectra  at  low  incident 
kinetic  energy.  The  extension  presented  here  is  necessary  when  the  orbitals 
of  various  bands  have  different  spatial  extent,  so  that  the  Auger  process  is 
controlled  not  only  by  the  density  of  states  but  also  by  the  relative  orbital 
length  scales. 

It  is  generally  true  that  the  two  electron  spectroscopic  methods  measuring 
convolutions  of  the  desired  quantities  are  at  a  disadvantage  when  compared 
to  their  one  electron  competitors,  since  some  of  the  details  might  be  lost  in 
the  deconvolution  process.  However,  if  interpreted  carefully  IN  can  be  a 
useful  complement  to  photoelectron  spectroscopy  due  to  its  extreme  surface 
sensitivity.? 

■\ 
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Most  other  spectra  have  a  low  energy  peak.(b>  The  Ni  3d  band  has  a  high 

density  of  states  and  a  short  spatial  extent,  while  the  s-p  band  has  a  low 

density  of  states  but  extends  further.  This  situation  should  cause  difficulties 
in  the  Hagstrum  model,  but  not  in  ours.  .A.  computation  will  establish  how 
serious  these  difficulties  are. 

The  parameters  used  in  the  calculation  are  given  in  Table  I.  Since  the 
constant  Cj^,C^,C^,C^  are  unknown  we  can  only  calculate  the  shape  of  the 
spectrum.  This  is  not  a  shortcoming  since  the  experimental  measurements  do 
not  give  the  absolute  intensities.  The  calculation  is  used  to  find  the 
density  of  states  that  best  fits  the  IN  data.  To  do  this  we  express  the 

density  of  states  as  a  sum  of  many  Gaussians  with  unknown  parameters  which 

are  varied  to  minimize  the  difference  between  the  computed  and  the  measured 
spectrum.  The  minimization  is  automated:  a  subroutine  generating  the  IN 
spectrum,  according  to  the  method  outlined  in  Section  II,  for  a  given  set  of 
Gaussian  parameters  is  used  in  a  Fletcher-Powell  minimization  program  which 
varies  the  parameters  until  the  best  least  square  fit  is  obtained. 

The  best  fit  of  the  Ni(lll)  IN  spectriun  is  shown  in  Figure  8.  The 
discrepancy  at  low  kinetic  energy  is  perhaps  due  to  the  presence  of  secondary 
electrons  in  the  experimental  spectrum.  In  Figure  9  we  show  the  best  density  of  states 
(full  line)  as  well  as  chat  obtained  with  the  Hagstrum  model.  In  Figure  10 
we  show  a  comparison  between  the  experimental  IN  spectrum  and  a  calculation 
using  Che  best  density  of  states  and  making  the  assumption  that  ion  neutralization 

9 

takes  place  only  at  the  ion  surface  distance  of  z=2.4A.  The  error  introduced 
by  this  assumption  is  large  enough  to  give  a  visibly  worse  fit  chan  the  full 
model. 

IV.  SUMMARY 

In  this  paper  we  have  developed  a  practical  method  for  the  analysis 
of  ion  neutralization  spectra.  This  was  stimulated  by  the  fact  chat  metastable 
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alone.  The  attractive  part  of  the  force  is  given  by  the  image  formula 

—20  -2 

F  =  (e-l)(£+l)  5.76x10  qz  (19) 

where  £  is  the  dielectric  constant  of  the  metal,  q  is  the  ion  charge  and  z  is  the 

O 

ion  surface  distance  in  A. 

II.  6.  Electron  escape  probability. 

Finally,  the  last  quantity  needed  for  the  calculation  of  the  Auger 
spectrum  is  the  escape  probability  which  we  use  a  semi-empirical 

formula  (Equation  (B.4)  of  Appendix  B)  proposed  by  Hagstrum,  with  the 
parameter  f  equal  to  2.2. 

III.  NUMERICAL  STUDIES  OF  THE  He''‘/Ni(lll)  SYSTEM. 

We  have  used  the  model  outlined  in  Section  II  to  analyze  the  IN  spectrum 
of  the  Ni(lll)  surface,  obtained  by  metastable  quenching  spectroscopy.^  The 
method  probes  the  surface  with  a  thermal  beam  of  neutral  He  and  excited 
(metastafale)  He*(2^S)  atoms.  The  metastable  is  resonantly  ionized  by  the  surface 
and  the  He"*"  ion  is  Auger  neutralized  when  it  collides  with  the  metal.  Since 
resonant  ionization  is  a  long  range  process  and  Auger  neutralization  is  a 
short  range  one  we  assume  chat  we  deal  with  a  two  step  process:  first 
resonant  ionization,  then  ion  neutralization.  In  this  case  the  resulting 
IN  spectrum  is  the  same  as  that  obtained  with  a  low  energy  (~0.05  eV)  ion  in 
beam. 

To  compare  Che  experimental  and  the  computed  spectra  we  must  remove  the 
secondary  electron  emissiai from  the  data.  In  spite  of  the  existence  of 
several  experimental  scudies^^  of  secondary  emission,  it  is  difficult  to 
perform  such  a  subtraction  accurately.  Because  the  precise  shape  and  magnitude 
of  the  low  energy  part  (£^  <  5  eV)  of  the  IN  spectrum  is  somewhat  uncertain,  it 
was  not  included  in  our  optimization  calculation. 

There  are  several  reasons  for  choosing  Ni(lll)  for  this  study:  (a)  its 
IN  spectrum  has  a  peculiar  form,  with  a  peak  at  high  electron  kinetic  energy. 
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exp(-Br)  ,  given  in  spherical  coordinates  with  the  origin  of  the 
coordinate  systeti  at  the  center  of  the  ion.  Since  the  ejected  Auger  electrons 
of  interest  in  IN  have  a  kinetic  energy  of  -13  eV  or  less,  their  Da  Broglie 
wavelength  is  larger  than  the  size  of  the  ionic  orbital  which  determines 
the  integration  range  in  Equation  (.IS)-  Therefore,  we  can  take  t:-o 

be  a  constant  C^.  The  screened  potential  is  taken  to  be  of  Yukawa  form 
W(ix^-x^l)  =  exp(-A.lx2-x^l ) )  /  )  •  We  have  also  used 

W(|X2-x^|)=  (|x2-x^1-a)  and  the  results  obtained  with  the  two  equations  are 
very  similar. 

The  calculation  of  the  matrix  element  is  outlined  in  Appendix  A.  The 
result  is 

(2Tr)^C  C  C  C 

-H(£.  ,£  .£  :2) - -  ^ - {exp(-6(2-Y)-(a  +X)  ] 

^  X(a^-X)6(6-a^-X)  ^ 

[(z+Y)  +  (B-a^-X)'^  +3'^]  (18) 

-expC-(a^+X)z]  [(B-a^-X)"^  +6~^]} 

This  ecuation  is  valid  for  a  >X>o,  8>(a  -X)  and  A>y  (see  .Apnendix  .A).  Here  Y 

V  ii 

is  a  cut-off  distance  which  is  the  ion's  closest  possible  approach  to  the 
metal  surface. 

II.  5.  The  trajector'/  of  the  incident  ion. 

.At  low  kinetic  energy  it  becomes  difficult  to  choose  a  proper  trajectory 
for  the  ion.  This  is  a  general  problem  for  all  curve  crossing  models. The 
system  is  simultaneously  on  curve  one  and  rwo,  with  a  given  transition  amplitude 
for  each  event.  Therefore  neither  energy  surface  is  providing  the  force  alone 
and  the  total  force  in  the  correct  equation  of  motion  for  the  ion  must  depend  on 
the  gradients  of  the  two  potential  energies  and  on  the  probability  amplitudes 
describing  the  occupation  of  the  two  electronic  states. In  spite  of  this  we 
use  here  an  equation  of  motion  corresponding  to  the  motion  on  the  ionic  surface 
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simplest  possible  model  in  which  the  orbital  lengths  of  various  bands  is 
qualitatively  incorporated. 

The  matrix  element  characterizing  the  Auger  neutralization  process  Is 
;z)  -  /  dr^  /  dr2  W(  j  r ^^-^2  1  )':^y (r j_)0,^(r2)  ,  ( 


12,13 


16) 


I  ^  1 

where  the  orbitals  were  defined  above  and  W(|r^-r2l)  is  a  screened 


Coulomb  interaction  between  two  surface  electrons. 


Since  electrons  are  indistinguishable  particles  we  should  subract  an 

]_3  *-►  *  *-► 

exchange  term  (which  is  obtained  by  replacing  with 

in  Equation  (16))  from  the  Coulomb  term  Included  in  Equation  (16).  Previous 

12 

surface  Auger  neutralization  work,  which  computed  these  matrix  elements,  has 
neglected  exchange.  The  reason  for  this  neglect  is  not  clear,  nor  is  it 
documented  that  the  term  is  negligible.  In  the  present  context,  where  the  main 
role  of  these  matrix  elements  is  to  take  into  account  the  fact  that  the  orbitals 
in  different  bands  have  different  spatial  extent,  the  inclusion  of  the  exchange 
term  does  not  seem  critical. 


In  order  to  model  the  qualitative  behavior  of  the  matrix  element 
H(e^,£^,£^;z)  we  use  the  following  simple  forms  for  the  orbitals.  The  orbitals 
and  describing  the  metal  atoms  are  taken  to  be 


Vg(r^)  =.  Cg  exp(-agZ3_) 


(17) 


where  the  constants  C,  and  have  different  values  whether  £,  is  in  band  A  or 

band  B.  In  the  case  of  Ni  the  two  bands  are  the  3d  band  and  the  s-p  one.  Since 

1-4,12 


there  are 


indications 


that  Auger  neutralization  takes  place  at  a 


relatively  large  distance  from  the  surface  (~2A) ,  and  that  at  those  distances 

14 

the  corrugation  of  the  electron  charge  density  is  very  low,  we  neglect  the 
change  in  ->3(r^)  along  the  surface.  The  ion  orbital  is  of  the  form 
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APPENDIX  B 

PROBABILITY  OF  ESCAPE  OF  AN  EXCITED  ELECTRON  FROM  A  METAL  SURFACE 


I 


The  probability for  an  excited  electron  of  energy  to  escape  from 


a  mecal  surface  is  given  by^’^ 

27100 

(P(e,)  =  /  /  (P(0  ,e  )sir6deilit 


(B  1) 


L 


where  maximum  value  of  9  for  which  escape  over  the  surface 

barrier  is  possible!  Here  measures  the  energy  of  the  ejected  Auger  electron 
relative  to  the  vacuum.  If  we  assume  that  the  angular  distribution  of 
electrons  is  isotropic  then  (P(9,£j^)  is  constant: 


(?(e,e.)  »  (477) 


-1 


(B  2) 


Substitution  of  Equation  (B  2)  into  Equation  (B  1)  gives 

•s 


l-tE  /(£,+€  for  t.  >  0 

0  k  o  k 


for  £,  <  0 
k 


(B  3) 


where  £  is  the  energy  of  the  bottom  of  the  conduction  band  with  respect  to 
o 

vacuum.  Calculations  of  the  total  electron  yield  using  the  Equation 
(B  3)  for(P(£^)  give  results  significantly  below  those  determined  experimentally. 
This  discrepancy  between  the  calculated  and  experimentally  determined  total 
.Auger  electron  yield  suggest  that  is  not  an  isotropic  function. 

To  model  the  anisotropy  in  the  Auger  electron  ejection  pattern,  the  following 
semi-empirical  expression  for  the  electronic  escape  probability  was  developed 
by  Hagstrum 


► 


The  numerical  value  of  Che  parameter  f  can  be  determined  by  fitting  the 
calculated  total  Auger  electron  yield  to  experimental  measurements.  For  helium 
ions  incident  on  a  clean  tungsten  surface  vich  kinetic  energy  of  40  eV , 

Hagscrum  gives  a  value  of  f  =  2.2.  Hagscrum^  finds  the  value  of(?(£j^)  for 
(£^)  >  4  eV  is  highly  insensitive  to  parametric  variation,  so  error  in  f  do 
not  lead  to  significant  errors  in  the  computed  IN  spectrum. 
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Table  1.  Physical  constants  used  for  the  He  /Ni(lll)  system. 


Helium  ionization  energy: 
ilatE 

Ca) 


(Isolated  atom) 


E.  (z_) 


m 


24.6  eV 

22.6  eV 


Exponential  decay  parameters: 


He  wave  function 

1  i 


8  =  3.779  A 


Ni(lll)  conduction  band  wave  functions : 


\s^^^  ”  ^4s® 


a,  -  0.35  A 
4s 


-1 


'^3d(2)  =  C3^e 


■^3d^' 


(c) 

Ni(lll)  conduction  band  energetics  : 
Fermi  level 
Lower  band  edge 


e„  =  -  5.2  eV 

f 

^LBE  ” 


(a)  Ionization  energy  at  the  position  of  most  probable  neutralization. 

26 

(b)  Decay  constant  of  atomic  nickel  wave  functions. 

3 

(c)  All  energies  measured  relative  to  the  vacuum  level. 


-26- 


LIST  OF  FIGURES 

1.  Energy  level  diagram  describing  (a)  the  initial  -and  (b)  the  final  states 
in  Auger  neutralization  of  the  incident  ion  by  a  metal  surface.  The  ion 
is  represented  by  the  half  filled  orbital.  The  metal  orbitals 

and  are  involved  in  Auger  neutralization.  The  ion  surface 
distance  is  z. 

Two  electronic  surfaces  describing  the  energy  of  the  ground  state  metal  and 
the  ion  (Me  +  A^) ,  and  the  energy  of  the  neutral  atom  A  and  the 

A 

metal  (Me  )^^j^having  an  electron  in  and  holes  in  and  The  orbital 

energies  are  defined  in  Figure  1. 

A  schematic  plot  of  the  neutral  curve  and  two  ionic  curves  corresponding 
to  two  values  of  ionic  curves  have 

different  asymptotic  energy  mismatches  A.  , («)  and  /L  "(»)  and  different 

cC^Ct  iCyCC 

crossing  posaitions  Z(a)  and  Z(a'). 

A  schematic  plot  of  the  family  of  ionic  curves  labelled  by  various  values 
of  a  =  which  can  participate  in  Auger  neutralization. 

A  schematic  description  of  the  bands  A  and  B  and  the  orbital  ranges 
Z^  and  Zg  discussed  in  the  text. 

A  schematic  plot  of  the  kinetic  energy  ranges  for  the  processes  BB,  .AB, 

AA  and  BA  defined  in  the  text. 

A  schematic  plot  of  the  IN  spectrum  for  two  situations.  Full  line:  The 

IN  spectrum  for  Z  >  Z  ,  £  >e„  and  P„>o  (Z  and  Z  are  defined  in  Figure  5. 
oAAooAdA 

£.  and  are  the  orbital  energies  in  band  A  and  B  and  and  c,  are  the 
A  o  D  A 

density  of  states  in  the  bands  A  and  B) .  Broken  line:  The  IN  spectrum 
for  2^>Zg,  and  The  AA,  AB,  BA  and  BB  processes  are  defined 

in  the  text.  The  energies  £p,£^,and  €2  defined  in  Figure  5.  I^  and 
I„  are  the  ionization  potentiaJs  of  the  atom  at  the  distances  Z.  and  . 


8.  Kinetic  energy  distributions  of  electrons  ejected  by  thermal  helium  ions 
colliding  with  a  clean  Ni(lll)  surface.  The  experimental  spectrum  was 
obtained  by  MQS.^  The  theoretical  spectrum  was  computed  by  using  the 
extended  model  described  in  Section  II.  The  density  of  states  was  varied 
to  get  the  best  fit  of  the  experimental  IN  spectrum. 

9.  The  density  of  states  of  Che  Ni(lll)  surface.  The  dotted  line  was  obtained 
by  using  the  extended  model  and  varying  the  density  of  states  to  produce 
the  best  fit  of  Che  experimental  IN  spectrum.  The  full  line  was  obtained 
by  using  Hagstrum  model. 

10.  The  IN  spectrum  of  the  Ni(lll)  surface,  computed  by  using  the  "best 
density  of  states"  shown  in  Figure  9  and  by  assuming  that  Che  neutralization 
of  Che  ion  takes  place  at  z  =  2.4  A. 

11.  The  coordinate  system  and  the  notation  used  for  the  evaluation  of  the 


matrix  elements. 
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Atom  -  surface  distance 


Electron  Kinetic  Energy 
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